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Summary. Current analyses of VLT/UVES quasar spectra disagree with the 
Keck/HIRES evidence for a varying fine-structure constant, a. To investigate this 
we introduce a simple method for calculating the minimum possible uncertainty on 
Aa/a for a given quasar absorber. For many absorbers in Chand et al. (2004) and 
for the single-absorber constraint of Levshakov et al. (2006) the quoted uncertain- 
ties are smaller than the minimum allowed by the UVES data. Failure of this basic 
consistency test prevents reliable comparison of the UVES and HIRES results. 



1 Introduction 

The 'many-multiplet' (MM) method is the most precise technique for con- 
straining cosmological changes in the fine-structure constant, a, from QSO 
absorption spectra [1]. We have previously described self-consistent MM ev- 
idence from 143 Keck/HIRES absorbers for a smaUer a over the redshift 
range 0.2 < Zabs < 4.2 at the fractional level Aa/a = (-0.57 ± 0.11) x 10~^ 
[2] . Clearly, independent analyses of spectra from different spectrographs are 
desirable to refute/confirm this. First attempts with VLT/UVES spectra, 
e.g. [3, 4], generally found null results with quoted uncertainties < 0.1 x 10~^. 
To investigate these claims, we introduce a simple method for calculating the 
minimum possible uncertainty on Aa/a from a given absorption system. 



2 A simple measure of the limiting precision on /\cx./a 

The velocity shift, Vj, of transition j due to a small relative variation in a, 
Aa/a <C 1, is determined by the g-coefficient for that transition. 



ia./aof-l] ^ ^«-2^^, (1) 



where ojqj & ujz,j are the rest-frequencies in the lab and at redshift z, ao is 
the lab value of a and az is the shifted value measured from an absorber at z. 
The MM method is the comparison of measured velocity shifts from several 
transitions (with different q-coefficients) to compute the best- fit Aa/a. The 
linear equation (1) implies that the error in Aa/a is determined only by 
the g-coefiicients (assumed to have negligible errors) and the total velocity 
uncertainty, integrated over the absorption profile, from each transition, (Jv,j' 
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5{Aala\^ = ^sJd, (2) 



for S ^ (^)~' and D ^ S^.^ {^y' [E, ^ (^)" 

This expression is just the solution to a least-squares fit of y = a + 6x to 
data {xi,yi), with errors only on the yi, where the intercept a is also allowed 
to vary; this mimics the real situation in fitting absorption lines where the 
absorption redshift and Aa/a must be determined simultaneously. 

The quantity a^j is commonly used in radial- velocity searches for extra- 
solar planets, e.g. [5], but is not normally useful in QSO absorption- line stud- 
ies. Most metal-line QSO absorption profiles display a complicated velocity 
structure and one usually focuses on the properties of individual velocity com- 
ponents, each of which is typically modelled by a Voigt profile. However, it is 
important to realize that Aa/a and its uncertainty are integrated quantities 
determined by the entire absorption profile. Thus, a^j should incorporate 
all the velocity-ccntroiding information available from a given profile shape. 
From a spectrum F[i) with 1-a error array crF[i), the minimum possible 
velocity uncertainty contributed by pixel i is given by [5] 

c x{i) [dF{i)/dX{i)] ■ ^ ' 

That is, a more precise velocity measurement is available from those pixels 
where the flux has a large gradient and/or a small uncertainty. This quantity 
can be used as an optimal weight, W{i) = [i7v(«)/c]^ , to derive the total 
velocity precision available from all pixels in a (portion of) spectrum, 

a^ = c[j:,W{i)]-'^\ (4) 

For each transition in an absorber, a^j is calculated from (3) & (4); the only 
requirements are the 1-a error spectrum and the multi-component Voigt pro- 
file fit to the transition's absorption profile. The latter allows the derivative 
in (3) to be calculated without the influence of noise. Once ayj has been 
calculated for all transitions, the uncertainty in Aa/a follows from (2). 

It is important to realize that the uncertainty calculated with the above 
method represents the absolute minimum possible 1-a error on Aa/a; the 
real error - as derived from a simultaneous x^-minimization of all parameters 
comprising the Voigt profile fits to all transitions - will always be larger than 
S{Aa/a)iim from (2). The main reason for this is that absorption systems 
usually have several velocity components which have different optical depths 
in different transitions. Equation (2) assumes that the velocity information 
integrated over all components in one transition can be combined with the 
same integrated quantity from another transition to yield an uncertainty on 
Aa/a. However, in a real determination of Aa/a, each velocity component 
(or group of components which define a sharp spectral feature) in one tran- 
sition is, effectively, compared with only the same component (or group) in 
another transition. Fig. 1 illustrates this point. 
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Fig. 1. Left: Simulation of two transitions in a multi-component absorber. Labelled 
are distinct stages of differential saturation in the two main spectral features. Right: 
The difference between the limiting precision, 5{Aa/ a)iiui, and the actual precision 
(as derived by x^-minimization analysis) varies strongly from stage to stage. In a 
single-component absorber the actual uncertainty tracks 5{Aa/ a)iii^. 

3 Application to existing constraints on ^ol/ol 

We have calculated 5{Aa/ a)\im for the absorbers from the three independent 
data-sets which constitute the strongest current constraints on Aa/a: (i) The 
143 absorbers in our Keck/HIRES sample [2]; (ii) The 23 Mg/Feii systems 
from UVES studied by [3]; (iii) The UVES exposures of the Zabs = 1-1508 
absorber towards HE 0515— 4414 studied by [4]. For sample (ii), the spectra 
were kindly provided by B. Aracil who confirmed that they have the same 
wavelength and flux arrays as those used in [3]. The main difference is the 
error arrays: our error array is generally a factor w 1.4 smaller than that used 
by [3] (H. Chand, B. Aracil, 2006, priv. comm.). We have confirmed this by 
digitizing the absorption profiles plotted in [3]. Thus, 5(zia/a)iim calculated 
using our spectra will be smaller than the value [3] would derive. The con- 
tinuum normalization will also be slightly different, but this has negligible 
effects on the analysis here. For sample (iii), the exposures were reduced us- 
ing a modified version of the UVES pipeline and their S/N matches very well 
those quoted by [4]. For samples (ii) & (iii) we use the Voigt profile models 
pubhshed in [3] & [4] to calculate 5{Aa/ a)\ira- For all samples, the atomic 
data for the transitions (including g-coefficients) were the same as used by 
the original authors. In practice, when applying (3) and (4) we sub-divide the 
absorption profile of each transition into 15kms~^ chunks to mitigate the ef- 
fects illustrated in Fig. 1. This provides a value of 5{Aa/ a)\n-a,k for each 
chunk k. The final value ol6{Aa/a)\\ra is simply {^j, l/[(5(Z\a/a)iim,fc]^}^^^; 
in all cases this is < 1.4 times the value obtained without sub-divisions. 

Figure 2 shows the la error on Aa/a quoted by the original authors ver- 
sus the limiting precision, 5{Aa/ a)\im. The main results are clear. Firstly, 
the Icr errors quoted for the HIRES sample in [2] always exceed 5(Z\a/a)iim, 
as expected if the former are robustly estimated. Secondly, for at least 11 
of their 23 absorbers, Chand et al. [3] quote errors which are smaller than 
5{Aa/ a)\ia\- Recall that since their error arrays are larger than ours, 11 out 
of 23 is a conservative estimate. Finally, the very small error quoted by [4] for 
HE 0515— 4414, 0.084 x 10~^, disagrees significantly with the limiting preci- 
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Fig. 2. Quoted errors vs. limiting precision for current samples. The Chand et 
al. and Levshakov et al. samples fail the basic requirement 5{Aa/a) > S{Aa/a)iim. 

sion of 0.14 X 10^^. Thus, the (supposedly) strong current UVES constraints 
on Aa I a fail a basic consistency test which not only challenges the precision 
claimed by [3] and [4] but which must bring into question the robustness and 
validity of their analysis and final Aa/a values. 

4 Conclusion 

We have introduced a very simple method of determining the limiting pre- 
cision on Aa/a obtainable from a set of transitions in a QSO absorption 
system, (5(Z\a/a)ii,n- Only the 1-cr error spectrum and the model absorp- 
tion profile (generally constructed from Voigt profiles) are required to cal- 
culate 5{Aa/ a)\ira. The method simply equates the total velocity informa- 
tion contained in the absorption profile of a transition with the expected 
error on Aa/a via that transition's sensitivity to a- variation. All absorption 
systems in our HIRES data-set have quoted uncertainties in Aa/a which 
exceed 5{Aa/a)\{m, as expected. However, the uncertainties on Aa/a cur- 
rently quoted by [3] and [4] from UVES spectra are, in many cases, smaller 
than 5{Aa/ a)\ira- Clearly, the UVES and HIRES results cannot be rehably 
compared until more rigorous analyses of the UVES spectra are completed. 
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